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Abstract 

A sizable rate of events, with several pairs of 6-quarks produced contemporar- 
ily by multiple parton interactions, may be expected at very high energies as 
a consequence of the large parton luminosities. The production rates are fur- 
ther enhanced in hadron-nucleus reactions, which may represent a convenient 
tool to study the phenomenon. We compare the different contributions to 
bbbb production, due to single and double parton scatterings, in collisions of 
protons with nuclei at the CERN-LHC. 
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I. INTRODUCTION 



The large rates of production of heavy quarks, expected at high energies, may lead to 
a sizable number of events, at the CERN-LHC, containing two or more pairs of b quarks, 
generated contemporarily by independent partonic collisions. An inclusive cross section 
of the order of 10 fib may in fact be foreseen for a double parton collision process, with 
two bb pairs produced in a pp interactions at 14TeV, while the cross section to produce 
two bb pairs by a single collision at the same cm. energy may be one order of magnitude 
smaller [1]. All production rates are significantly enhanced in proton- nucleus collisions, 
which may offer considerable advantages for studying multiparton collision processes [2]. 
Given the large rates expected, the production of multiple pairs of 6-quarks should hence 
represent a convenient process to study multiparton interactions in pA collisions at the 
LHC. On the other hand the mechanism of heavy quarks production is not yet understood 
satisfactorily also in the simplest case of nucleon-nucleon collisions, the effects of higher order 
corrections in a s being still a matter of debate. A comprehensive description, of the much 
more structured process of heavy quarks production in hadron-nucleus interactions, might 
hence be approached after gaining a deeper understanding of the short scale parton-level 
dynamics of the process. On the other hand a significant feature of higher order corrections in 
a s is that, for a limited set of physical observables, the whole effect of higher orders reduces 
to an approximate rescaling of a lowest order calculation of heavy quarks production in 
perturbative QCD [3]. Some of the features of the process are therefore effectively described 
by the simplest parton level dynamics at the lowest order in a s , which lets one speculate that 
a similar property might hold also for a much more complex process, as bbbb production. 

Taking this optimistic point of view we will attempt, in the present note, to obtain indi- 
cations on some properties of bbbb production in pA interactions at the LHC, by considering 
the contributions of the two different interaction mechanisms, the connected 2 — > 4 and 
the disconnected (2 — > 2) 2 parton processes, which will be effectively described by the low- 
est order diagrams in perturbative QCD, keeping higher order corrections into account by 
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a simple overall rescaling. Obviously, following this philosophy, all considerations will be 
limited to the class of physical observables falling in the category mentioned above. 

II. BBBB PRODUCTION IN PROTON-NUCLEUS COLLISIONS 

Quite in general [4] [5], with the only assumption of factorization of the hard component 
of the interaction, the expression of the double parton scattering cross section to produce 
two bb pairs is given by 

ct (p,a)( & M6) = \ Y1 J ^p{xi,xf,Si j )a(xi,x[)a(x j ,x' j )Ti yN) A){x' j ,x'^Si j )dx i dx[dx (1) 

ij 

where the index N refers to nucleon and A to a nucleus, while the indices i,j to the different 
kinds of partons that annihilate to produce a bb pair and the factor 1/2 is a consequence of the 
symmetry of the expression for exchanging % and j. The interaction region of a hard process 
is very small as compared to the hadron scale, so the two interactions are well localized in 
transverse space within the two overlapping hadrons. The rate of events where two hard 
collisions take place contemporarily, in a given inelastic hadron-hadron interaction, depends 
therefore on the typical transverse distance between the partons of the pairs undergoing the 
multiple collision processes. A main reason of interest is hence that double parton scatterings 
may provide information on the typical transverse separation between pairs of partons in the 
hadron structure. Indeed the non-perturbative input of a double parton collision is the two- 
body parton distribution function T(xi, x 2 , Si j2 ), which depends not only on the fractional 
momenta x± t 2, but also on the relative distance in transverse space s± t 2, besides (although 
not written explicitly to simplify the notation) on the scales of the two interactions and on 
the different kinds of partons involved. As a consequence the double parton scattering cross 
section is characterized by a linear dependence on dimensional scale factors, which have a 
direct relation with the typical transverse distances between the partons of the various pairs, 
contributing to the double scattering process under consideration. 

The cross section is simplest when the target is a nucleon and partons are not corre- 
lated in fractional momenta, which may be a sensible approximation in the limit of small 
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x. In such a case the two-body parton distribution may be factorized as T p (xi,Xj] Sij) = 
G(xi)G(xj)F(sij), where G(x) are the usual one-body parton distributions and F(s) a func- 
tion normalized to 1 and representing the parton pair density in transverse space. The 
inclusive cross section to produce two pairs of b quarks is then written as [6] 



^£(66; 66) = -]T 6^(66)^(66) 



(2) 



where (Jj(66) represents the inclusive cross sections to produce a bb pair in a hadronic collision, 
with the index % labelling a definite parton process. The factors B*- 7 have dimension an inverse 
cross section and result from integrating the products of the two-body parton distributions 
in transverse space. In this simplified case the factors iJ provide a direct measure of 
the different average transverse distances between different pairs of partons in the hadron 
structure [6] [7]. 

The cross section has a more elaborate structure in the case of a nuclear target. The 
most suitable conditions to study the phenomenon are those where the nuclear distributions 
are additive in the nucleon parton distributions. In such a case one may express the nuclear 
parton pair density, 1^4(2^,2^; s^), as the sum of two well defined contributions, where the 
two partons are originated by either one or by two different parent nucleons: 



r^rr,, Xj] Sij) — T^(x i , Xj] s^j) + 1^4(2^, Xj] s 



'■ij, 



(3) 



and correspondingly — <x^|i + <r^| 2 . The two terms ^a\i,2 are related to the nuclear 
nucleon's density. Introducing the transverse parton coordinates B ± where B is the 
impact parameter of the hadron-nucleus collision, one may write 



TaM, x' p Sy )| 1>2 = / d 2 B lA (x[, x>; B+ S f,B- S f) 



1,2 



where 7a 1 1,2 are given by 

lAfa,^,B+ S f,B- S f) 



i = r N (x' i ,x' j ; Sij )T(B) 



^G^G^x'^B+'f^B-f 



(4) 



(5) 



with T(B) is the nuclear thickness function, normalized to the atomic mass number A and 
Gn nuclear parton distributions divided by the atomic mass number. 

In the simplest additive case, the first term in Eq.(3) obviously gives a simple rescaling 
of the double parton distribution of an isolated nucleon 

T A (x' t ,x' r , = rV^,4; s tj ) J d 2 BT(B) (6) 

and the resulting contribution to the double parton scattering cross section is the same as 
in a nucleon-nucleon interaction, apart from an enhancement factor due to the nuclear flux, 
which is given by the value of the atomic mass number A: 



Aa D N . (7) 



The cr®\2 term has more structure. In this case the integration on the relative transverse 
distance between the partons of the interacting pairs, s^-, involves both the projectile and 
two different target nucleons: 

J ds ij r p (x i ,x j ;s ij )T(B+ S f)T(B- S f). (8) 

As one may notice in this case two very different scales are compared, the hadron radius 
r p and the nuclear radius Ra- A usual approximation in pA interactions is to consider the 
limit r p <C Ra, where one may use the approximation 

T(S±^)~T(S), (9) 

which allows one to decouple the integrations on s^- and on B. One hence obtains: 

- J G p (xi,Xj)a(xi,x' i )a(xj,x' j )GN(x' i )GN(x'j)dx i dx' i dxjdx' j J d 2 BT 2 (B), (10) 
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where 



Gp{xi,Xj) J d SijY p^Xi, x j , Sy) • (l-Q 

Remarkably the two terms a® and a® have very different properties. In fact the 
correct dimensionality of a a is provided by transverse scale factors related to the nucleon 



scale, cfr. Eq.s(2),(7). The analogous dimensional factor in a% is provided by the nuclear 
thickness function, which is at the second power, being two the target nucleons involved in 
the interaction. 

As pointed out in ref . [2] , while on general grounds a® p A ^ depends both on the longitu- 
dinal and transverse parton correlations, the a® 2 term depends solely on the longitudinal 
momentum fractions X ^ j X j SO that, when the term is isolated, one has the capability 
of measuring the longitudinal and, a fortiori, also the transverse parton correlations of the 
hadron structure in a model independent way. 

Although the two contributions may be defined also in a more general case, the separation 
of the cross section in the two terms a® and a® is most useful in the regime of additivity 
of the nuclear structure functions, which may not be a bad approximation for a sizable part 
of the kinelatical regime of bottom quarks production at the LHC. In the case of a central 
calorimeter with the acceptance of the ALICE detector (|^| < 0.9), the average value of 
momentum fraction of the initial state partons, in a pp — > bbbb process, is (x) pa 6 x 10~ 3 . 
By introducing a cut in the transverse momenta of the b quarks of 5 GeV one obtains 
(x) ^ 1CT 2 , while a cut of 20 GeV in p t) within the same pseudorapidity range, gives 
(x) pa 2 x 10~ 2 . If considering a more forward detector, as LHCb (1.8 < rj < 4.9), the 
average value of momentum fraction is (x) pa 5 x 10~ 2 . Deviations from additivity at low 
x are less than 10% for x > 2 x 10~ 2 [8] and, although increasing with the atomic mass 
number, non additive corrections are at most a 20% effect, on the kinematical regime above. 

As already mentioned, heavy quarks production is characterized by a non trivial dy- 
namics, in such a way that also next-to-leading corrections to the lowest order term in 
perturbative QCD are not sufficient for an exhaustive description of the inclusive spectra, 
which most likely need an infinite resummation to be evaluated. By comparing the results 
of different approaches to heavy quarks production (as NLO QCD and fc r factorization) one 
nevertheless finds that, in a few cases, the whole effect of higher orders is to a large extent 
just an approximate rescaling of the results obtained by a lowest order evaluation in per- 
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turbative QCD. When limiting the discussion to an accordingly restricted set of physical 
observables, the whole effect of higher order corrections to heavy quark production is hence 
summarized by a single number, the value of the so-called K-factor. By evaluating, with the 
/ct-factorization approach, rapidity and pseudorapidity distributions of bb production, in pp 
collisions at the center-of-mass energy of = 5.5 TeV, within \r)\ < 0.9 and 1.8 < r] < 4.9, 
with different cuts in the transverse momenta of the produced b quarks, one finds a result 
not incompatible with a lowest order calculation in perturbative QCD rescaled by factor 
K = 5.5. 

To discuss the production of bbbb in pA collisions, while remaining in a kinematical 
regime where non additive corrections to the nuclear structure functions are not a major 
effect, we have limited all considerations to physical observables, where higher orders may 
be taken into account by the simple rescaling of the lowest order calculation. For bbbb 
production, where only results of three level calculations are up to now available, we have 
further assumed that the value of the if-factor in the 2^4 parton process is the same as in 
the 2^2 process. We have hence evaluated the various contributions to the cross section in 
the case of a central calorimeter with the acceptance of the ALICE detector, \r)\ < 0.9, and 
of a more forward detector as LHCb, 1.8 < rj < 4.9. In both cases the typical values of x 
are rather small, so the factorization of the cross section in Eq.(2) may not be unreasonable. 
Since the process is dominated by gluon fusion, the expression of the cross section in Eq. (2) 
may be limited to a single term only. For the corresponding dimensional scale factor we 
have used the value reported by the CDF measurement of double parton collisions [9] [10]. 
To obtain the cross section at the lowest order in pQCD we have used the MRS99 parton 
distributions [11], with factorization and renormalization scale equal to the transverse mass 
of the 6-quark. The cross sections of the 2^4 processes have been evaluated by generating 
the matrix elements of the partonic amplitudes with MadGraph [12] and HELAS [13]. For 
the mass of the bottom quark we have used the value m& = 4.6 GeV. The multi-dimensional 
integrations have been performed by VEGAS [14] and all lowest order pQCD cross section 
have been finally multiplied by a factor K = 5.5. 
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III. RESULTS 



A major result of the present analysis is that the effects induced by the presence of the 
nucleonic degrees of freedom, in double parton scatterings with a nuclear target, cannot 
be reduced to the simple shadowing corrections of the nuclear parton structure functions, 
whose effect is to decrease the cross section as a function of A. In the case of double parton 
collisions, the main effect of the nuclear structure is represented by the presence of the a^\ 2 
term in the cross section, which scales with a different power of A as compared to the single 
scattering contribution, producing an additive correction to the cross section. 

To emphasize the resulting "anomalous" dependence of the double parton scattering 
cross section, as a function of the atomic mass number, we have plotted in fig.l the ratio 
a 2 7( a F + aS )i as a function of A. The ratio represents the contribution to the cross section 
of the processes where two different nuclear target nucleons are involved in the interaction, 
scaled to the contribution where only a single target nucleon is involved. The dependence 
on the atomic mass number of the latter terms of the cross section, namely the single (<j s ) 
and the double (cxf) parton scattering terms against a single nucleon in the nucleus, is the 
same of all hard processes usually considered, where nuclear effects may be wholly absorbed 
in the shadowing corrections to the nuclear structure functions. The contribution to the 
cross section of the o\f term is, on the contrary, "anomalous" , involving two different target 
nucleons in the interaction process. The ratio above hence represents the relative weights of 
the "anomalous" to the "usual" contributions to the double parton scattering cross section on 
a nuclear target. The plots in fig.l refer to the case of a forward calorimeter (1.8 < rj < 4.9) 
and of a central calorimeter < .9) with different cuts on the transverse momenta of the 
produced 6-quarks (p cut = 0, 5, 10 GeV/c). 

The different contributions to the cross section for bbbb production, due to interactions 
with a single or with two different target nucleons, are shown in fig. 2 in the case of a central 
calorimeter. The left figure shows the cross section as a function of the atomic mass number 
of the one-nucleon (dashed line) and of the two nucleons (continuous line) contributions 
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to the cross section. The dotted line is the sum of the two terms. The right figure shows 
the two contributions to the rapidity distribution of a 5-quark produced in an event with 
two bb pairs: one-nucleon (dashed histograms) and two-nucleon contributions (continuous 
histograms) in the case of a heavy, Au, (higher histograms) and of a light nucleus, 0, (lower 
histograms). The dependence on rj of the two histograms is essentially the same, showing 
that the effect of the single parton scattering term (the 2^4 parton process) is negligible 
in this kinematical regime. 

The A-dependence of the two different contributions, as a function of A, are shown in 
fig. 3 in the case of a central calorimeter, after applying a cut of 5 GeV/c (left figure) and 
of 10 GeV/c (right figure) in the transverse momenta of each produced 6-quark. Dashed, 
continuous and dotted lines have the same meaning as in fig. 2. 

The case of a forward calorimeter, 1.8 < rj < 4.9, is shown in fig. 4, where the dashed, 
continuous and dotted lines have the same meaning as before. 

Summarizing the large size of the cross section of bbbb production in hadron-nucleus 
collisions at the LHC (the values are of the order of one hundreds of /ib,) suggests that the 
production of multiple pairs of 6-quarks is fairly typical at high energies, hence representing 
a convenient channel to study multiple parton interactions. A rather direct feature, which 
is a simplest prediction and then a test of the interaction mechanism described above, is 
the "anomalous" dependence on A. The effects induced by the presence of the nucleonic 
degrees of freedom in the nuclear structure are in fact not limited, in this case, to the 
shadowing corrections to the nuclear structure functions usually considered, which cause a 
limited decrease (not larger than 20%, in the kinematical regime considered here) of the 
cross section for a hard interaction in hadron-nucleus collisions. When considering double 
parton scatterings, all nuclear effects can be exhausted in the shadowing corrections to 
the nuclear structure functions only in the <Ja\i term. The dominant effect of the nuclear 
structure is however due to the presence of the term in the cross section, which scales 
with a different power of A as compared to single scattering term, giving rise to a sizably 
larger correction, with opposite sign as compared to the shadowing correction, namely to an 
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increase of the cross section, which may become larger than 100% for a heavy nucleus. 
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FIG. 1. Relative weights of the terms with "anomalous" and "usual" A-dependence in the double 
scattering cross section for bbbb production. 
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FIG. 2. Different contributions to the cross section for bbbb production, in a central calorimeter. Left 
figure: cross section as a function of A of the one-nucleon (dashed line) and of the two nucleons (contin- 
uous line) contributions. The dotted line is the sum of the two terms. Right figure: differential rapidity 
distributions of a 6-quark produced in an event with two bb pairs. One-nucleon (dashed histograms) and 
two-nucleon contributions (continuous histograms) in the case of a heavy (higher histograms) and of a light 
nucleus (lower histograms). 
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FIG. 3. Different contributions to the cross section for bbbb production in a central calorimeter as a 
function of A, after applying a cut of 5 GeV (left figure) and of 10 GeV (right figure) in the transverse 
momenta of each produced &-quark: one-nucleon (dashed line), two nucleons (continuous line) processes and 
total (dotted line). 
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FIG. 4. Different contributions to the cross section for bbbb production, in a forward calorimeter. Left 
figure: cross section as a function of A of the one-nucleon (dashed line) and of the two nucleons (contin- 
uous line) contributions. The dotted line is the sum of the two terms. Right figure: differential rapidity 
distributions of a 6-quark produced in an event with two bb pairs. One-nucleon (dashed histograms) and 
two-nucleon contributions (continuous histograms) in the case of a heavy (higher histograms) and of a light 
nucleus (lower histograms). 
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